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Tetrahedral correlations in ^"Zr and ^^Zr 
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Axial octupole and tetrahedral correlations in ^''Zr and ®^Zr have been investigated using the 
generator coordinate method, applied to a basis generated by Skyrme HF+BCS calculations. We 
focus on the possible presence of states with tetrahedral symmetry and their stability with respect 
to octupole vibrations. We show that pairing significantly reduces the stability of the tetrahedral 
configuration and that a shallow mean-field tetrahedral minimum coexists with an axial octupole 
minimum. The contributions to the correlation energies coming from the tetrahedral degree of 
freedom and octupole axial deformation are discussed. 

PACS numbers: 27.60.+J, 27.50.+e, 21.10.Re 



I. INTRODUCTION 

The possible tetrahedral instability of finite niany- 
fermion systems has been predicted in Ref. fll| . Recently, 
using the similar arguments, it has been conjectured that 
the intrinsic ground- or a low- lying isomeric state of some 
nuclei may possess a tetrahedral deformation 0, Q • This 
kind of symmetry is rather common in molecules and 
metallic clusters [j] where the mutual geometric arrange- 
ment of the ions determines the shape of the system. In 
atomic nuclei, the situation is more complicated and only 
sufficiently large shell effects can generate deformed con- 
figurations. Although there is convincing theoretical ev- 
idence for the presence of tetrahedral states in nuclear 
systems, they have not been confirmed experimentally. 

Tetrahedral shapes are invariant with respect to the 
transformations of the point group T^ ■ In nuclei, they 
are realized at first-order through non-axial octupole de- 
formations of the nuclear density corresponding to the in- 
trinsic octupole moment Q32 oc r'^{Y-ii + ^3-2)- The key 
argument in favor of the stability of tetrahedral shapes 
is a direct consequence of the theory of point groups. 
The group T^ possesses 2 two- and one four-dimensional 
irrep. The unusual (in the context of nuclear structure) 
family of four-fold degenerate levels lead to a bunching of 
single-particle states resulting in rather large shell gaps 
and increased stability of specific configurations [2, S Q 

Among the predicted "tetrahedral magic" numbers are 
the Z = 40 proton and A^ = 40, 56 — 58 neutron num- 
bers Q. The Zr isotopes with 40 and 58 neutrons are 
thus predicted to be doubly-magic with respect to this 
symmetry. However these nuclei may also exhibit other 
types of octupole deformations: the coupling between the 
neutron d^/2 and ft-n /2 orbitals and the proton ^3 /2 and 
(79/2 orbitals can lead to both axial and non-axial oc- 
tupole correlations [6]. Moreover, octupole deformations 
are in competition with the quadrupole mode which may 
obscure the signature of the tetrahedral symmetry. Al- 



though experimental data is rather poor for ^'^Zr, the 
rotational band built on the ground state suggests the 
presence of a large quadrupole deformation ''/| . The evo- 
lution of shapes for the Zr isotopes above the closed shell 
nucleus ^°Zr is rather complex. The low energy spectrum 
of ^^Zr exhibits a pattern typical of a spherical nucleus, 
while 100-1042]. i^ave very large quadrupole deformations 
in their ground state together with co-existing low-lying 
oblate and spherical minima 18|. ^^Zr lies at the bor- 
der between these two regions and exhibits a transitional 
character, as demonstrated both experimentally [9] and 
theoretically @, ^Q, [III . In particular, the large experi- 
mental EG transition between the first excited O"*" state 
and the ground state suggests a strong mixing between 
co-existing shapes. This nucleus appears thus to be par- 
ticularly rich in terms of the different deformation modes 
that are in competition: spherical, oblate, prolate, tetra- 
hedral and axial-octupole shapes. 



The purpose of the present article is to analyze for 
the first time the role of tetrahedral configurations in the 
collective excitations of these nuclei. It is organized as 
follows. First, we investigate the competition between 
axial and tetrahedral octupole shapes using the Skyrme 
Hartree-Fock + BCS (HFBCS) approach to probe the po- 
tential energy landscape in the octupole directions. We 
then explore dynamical effects beyond mean-field, parity 
restoration and quantum fluctuations, in order to deter- 
mine whether states with large tetrahedral correlations 
are present at low excitation energy and thus whether 
there is a possibility to identify them experimentally. 
This analysis, performed using the generator coordinate 
method JJJ(GCM), aims at paving the way to a more 
comprehensive quantum treatment of all quadrupole and 
octupole degrees of freedom of the nuclear surface at the 
same time. The method has already been applied to the 
study of the excitation modes of the superdeformed Hg 
and Pb isotopes p^ . 



II. MEAN-FIELD CALCULATIONS 



in previous applications, we set po = O.lBfm 



In order to investigate the variation of the energy of 
nuclei as a function of several shape degrees of free- 
dom, the Hartree-Fock (HF) equations have been solved 
by discretization on a 3-dimensional mesh in coordinate 
space [lj|. Contrary to the usual way of solving mean- 
field equations in a truncated oscillator basis, this tech- 
nique has the advantage that it allows to describe nuclear 
configurations with any kind of shape with the same, 
high, numerical accuracy. In particular, this method has 
been used to describe nuclei from their spherical ground 
state up to fission with an accuracy of a few tenth of keV 
on energy differences [IJ, [l^ • 

Since to study a nucleus as a function of several shape 
degrees of freedom is computationally a heavy task, we 
have imposed a symmetry condition which simplifies the 
problem while retaining the main degrees of freedom 
relevant for a study of tetrahedral shapes. The mean- 
field density is required to be symmetric with respect to 
two mutually perpendicular planes. Such a constraint 
reduces the complexity of the problem four times but 
still allows to study the variation of the nuclear en- 
ergy as a function of all octupole degrees of freedom, 
although odd-?7i and even-m modes cannot be treated 
simultaneously [13|. The coupling between multipole 
moments with even and odd m-values is beyond the 
scope of our study, which is focused on the tetrahedral 
mode. Note that deformations corresponding to all mul- 
tipole moments with even-m values, in particular triax- 
ial quadrupole deformations are automatically taken into 
account by our method. 

We have performed a full set of calculations with two 
parametrizations of the Skyrme interaction : SIII, which 
was used in previous studies of the Zr isotopes with the 
same method [HI [3, and SLy4 [l3|. The results that 
we present in details correspond to the SIII force. Both 
forces predict very similar behavior of the energy as a 
function of octupole degrees of freedom around the spher- 
ical configuration. However their predictions for the en- 
ergy as a function of the axial quadrupole moment are 
qualitatively different. The main difference between both 
forces is that, at the mean-field level of approximation, 
the ^''Zr ground state obtained with SLy4 is spherical 
and the deformed configuration is excited by 3-4 MeV, 
whereas for SIII the two configurations are almost de- 
generate. Also the spherical configuration of ^^Zr forms 
a local minimum with respect to the quadrupole deforma- 
tion for SLy4 force, whereas it turns out to be unstable 
for SIII. 

The pairing interaction has been treated in the BCS 
approximation including the Lipkin-Nogami (LN) correc- 
tion [l3|. A zero-range density-dependent pairing inter- 
action has been used: 

V,a^r = \g^{l - P.)S{r - r') (^1 - ^) , (1) 

where i = n,p for neutrons and protons, respectively. As 



TABLE I: The pairing strengths Qi used for calculations and 
the averaged pairing gap in the ground state. 


SIII 


5p (MeV fm^) g 


„ (MeV fm^ 


Ap (MeV) 


A„ (MeV) 


»OZr 
9«Zr 


1100 
1050 


1300 

725 


1.048 
0.912 


1.415 
0.655 



The strengths of the pairing force for both nuclei are 
listed in table I together with pairing gaps A^ calculated 
as an average over the single-particle states within a 5 
MeV window around the Fermi level. These values repro- 
duce the "experimental" pairing gaps A, extracted from 
the odd-even mass staggering using a three-point filter 
from Ref. [l^. Since there are large uncertainties in the 
experimental gaps, especially around ^"Zr, we have also 
performed calculations with reduced pairing strengths, in 
particular with the values used in Ref. ^20] which give 
gaps twice smaller than the "experimental" ones. 
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FIG. 1: Variation of the energy as a function of the 
quadrupole Q20 moment calculated in the HFBCS approach. 

The HFBCS+LN calculations presented in Fig.[T]show 
the variation of the energy as function of the axial 
quadrupole moment. For *°Zr, the spherical and de- 
formed configurations are almost degenerate, while for 
^*Zr, the SIII interaction does not predict a spherical 
minimum. However calculations with the SLy4 force give 
a shallow spherical minimum with a depth of 200 keV. 
These differences between the two parametrizations are 
due to the transitional character of ^^Zr which makes the 
study of its quadrupole properties very sensitive to the 
details of the effective interaction. 

In order to probe the octupole susceptibility of the 
spherical configurations, we have performed two sets 
of calculations along the axial octupole path {Q32 was 
kept equal to zero) and the tetrahedral path (Q30 was 
kept equal to zero). In both cases, the quadrupole mo- 
ment was constrained to zero. The results are shown 



in Fig. [21 The octupole moments have been expressed 
through the dimensionless deformation parameters us- 
ing the relation: /^so = (Q3o)/Co, /332 = {Q32)/C2, 
3 

where Co = --—A^r\,C2 = Cq/V2 with tq = 1.2fm. 

Ait 
For each octupole mode, we have performed calculations 

with pairing switched off and on. In the HF approxima- 
tion, the spherical configuration is always unstable with 
respect to the octupole modes. The energy gained by 
octupole deformations is around 1 MeV for Y30 and is 
still larger for Y^2, around 2 MeV for ^°Zi and 3 MeV 
for ®^Zr. This result is consistent with calculations per- 
formed within the macroscopic-microscopic method with 
a Woods-Saxon potential Q . When pairing correlations 
are taken into account, the effects of octupole correla- 
tions are strongly reduced. For ^^Zr the energy gain 
does not exceed 1 MeV and for ^°Zr the octupole minima 
even disappear completely. Nevertheless the susceptibil- 
ity towards the ¥32 mode remains slightly larger than 
for the axial octupole mode. Note that octupole minima 
in ^*Zr are saddle-points. Very similar results concern- 
ing the tetrahedral instability of the spherical configura- 
tion of ^°Zr have been found in Skyrme HFBCS ^ and 
HFB [20, [22] calculations. The mechanism by which pair- 
ing correlations makes the tetrahedral minimum weaker 
is rather simple. In these nuclei, the single-particle gaps 
are not large in the spherical configuration and arc in- 
creased almost twice by tetrahedral deformations. As 
a result, shell effects favor the creation of a tetrahedral 
minimum. On the other hand, the pairing interaction 
is the largest for the spherical configuration and is very 
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FIG. 2: Energy as function of the octupole moments. The 
quadrupole deformation has been constrained to zero. The 
solid and dotted lines denotes the HFBCS results and pure 
HF results, respectively. Lines with solid and empty circles 
denote the energy as a function of the Q32 and Q30 moment, 
respectively. 



small in the tetrahedral minimum, resulting in a flat po- 
tential energy curve. 

This result shows that a static tetrahedral configura- 
tion may appear only as a result of a delicate balance 
between pairing and shell effects, which makes its predic- 
tion strongly dependent on the effective interaction that 
is used. Moreover should there be a static tetrahedral 
minimum, it is likely to be shallow and could be destroyed 
by dynamical correlations. One must also emphasize that 
our mean-field calculations were performed around the 
spherical configuration, by constraining the quadrupole 
moment to zero. The octupole minima that are obtained 
could therefore also be unstable against quadrupole de- 
formations. Yet, dynamical factors may also play in favor 
of octupole minima. In particular, when octupole corre- 
lations are included, the mean-field wave function breaks 
parity and parity projection brings a correlation energy 
which favors octupole correlation, especially for negative 
parity states. We will study these effects in the next sec- 
tion to determine the stability of a tetrahedral state with 
respect to correlations beyond a mean-field approach. 



III. OCTUPOLE DYNAMICS BEYOND A 
MEAN-FIELD APPROACH 

The generator coordinate method (GCM) allows at 
the same time to study the collective dynamics of a nu- 
cleus with respect to a collective variable and to restore 
the symmetries broken in a pure mean-field approach 
(see [131 ^-iid references therein). It is perfectly suited 
to our goal which is to determine a possible infiuence of 
the tetrahedral mode on low-lying collective states. 

The information given by a mean-field potential en- 
ergy surface is of a purely static nature and one must go 
beyond a mean-field approach to study its stability with 
respect to vibrations in the potential well. In cases where 
the energy surface presents several co-existing minima 
separated by low barriers, the very concept of a mean- 
field state with a specific deformation may even break 
down. Correlations beyond mean-field lead then to states 
which represent a mixture of configurations with differ- 
ent deformations. All these effects are taken into account 
in the GCM approach. 

We will limit here the collective space to octupole 
modes against which ^"Zr and ^*Zr are very soft. It is 
clear that a full study of these nuclei would require to 
include also quadrupole deformations and to mix both 
octupole modes. Such multi-dimensional GCM calcula- 
tions are feasible but represent a heavy task that is be- 
yond the scope of this exploratory article. The purpose 
of the present GCM calculations is limited to find out 
whether there is a chance to have a clear signature of the 
tetrahedral mode in the low energy spectrum. 

We have taken either Q30 or Q32 as generator coordi- 
nate. The wave- functions generated by the constrained 
mean-field calculations have first been projected on par- 



tide numbers and parity: 



i?(iV,Z,/33^)j 



(0(/33p)|i^P( 



±,Af,Z' 



\HM) 



{<f>{(33,m±,N,zMP3^)) 



(2) 



where \4>{fi3fi)) are HFBCS wave functions generated 
with the constraint ((/'(/33^)|(33mI0(/53m)) = C'm/^Sm- The 
operator P{±^n,z) is the product of operators projecting 
on TT = ±1 parity and on N neutrons and Z protons. The 
parity-projected energies are shown in the upper part 
of Fig. [31 As usual after parity restoration [13] the en- 
ergy minima for positive parity states are shifted towards 
smaher octupole deformations compared to the HFBCS 
minima (see Figl2]) while the negative parity states have 
larger deformations. For both nuclei the energy minima 
for positive parity correspond to very similar /S^q and 
/332 values. For the negative parity curve, /?32 is system- 
atically larger than P^q in the minimum. Qualitatively 
similar results have been obtained with the SLy4 Skyrme 
parametrization. 

The GCM allows to study the stability of the config- 
urations, corresponding to the minima of these energy 
curves, with respect to large amplitude vibrations. A col- 
lective wave function is constructed by mixing the mean- 
field states corresponding to different values of the oc- 
tupole moment, after their projection on particle number 
and parity: 



I*>=//(/33m)A 



{±,N,Z) 



\HP3^))df33^ 



(3) 



Er 



E{N,Z,spher.)-E^ 



(4) 



where E(N, Z, spher.) is the energy of the particle num- 
ber projected spherical configuration obtained in the HF- 
BCS approach, and E'^ is the lowest positive-parity en- 
ergy obtained in the GCM. We calculated also a dy- 
namical deformation associated with the lowest collective 
eigenstates for both parities defined by: 



/3; 



'3m 



03^ 



(5) 
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The coefficients ](Ji-ip^ are determined by minimization 
of the total energy of the collective wave function |^). 
The same effective interactions as in the mean-field cal- 
culations are used. In practice, the integral is replaced 
by a discrete summation over /33^, with a number of 
points large enough to obtain results independent of the 
discretization [ij|. The discretized Hill- Wheeler (HW) 
equation was solved separately for each collective coor- 
dinate Qso and Q32. The collective wave functions (re- 
lated to /(/33/i) by an integral transformation) are plot- 
ted in Fig. 3. One can see that these wave functions are 
spread around the minima of the projected mean-field 
energy curves, with a shape typical of a vibration in a 
1-dimensional energy well. 

The energy gain due to correlations beyond the mean- 
field is shown in table II. The correlation energy is defined 
by: 
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FIG. 3: The parity projected mean-field energy as function 
of the octupole moments is shown in the upper figures. The 
solid and dotted lines denote the positive and negative parity 
results, respectively. Lines with open and filled circles de- 
note the energy as a function of the Q30 and Q32 moment, 
respectively. In the lower figure, the square of collective wave 
functions are shown. The line styles correspond to those in 
the upper figure. 



TABLE II: Results of the GCM calculations. E^xc denotes the 
excitation energy of the first negative-parity collective state 
with respect to the first positive-parity state. Ecorr is the 
correlation energy as defined in (|4]) and /^a^j refers to the dy- 
namical deformation of ([5]). 



SIII 


£e.c (MeV) 


E.orr (MeV) 


TV 


/330 


A2 


»OZr 


0.0 


1.499 


+ 1 


0.0 


0.12 




2.22 


- 


-1 


0.0 


0.25 




0.0 


1.459 


+ 1 


0.13 


0.0 




2.52 


- 


-1 


0.23 


0.0 


8«Zr 


0.0 


1.485 


+ 1 


0.0 


0.14 




0.784 


- 


-1 


0.0 


0.21 




0.0 


1.387 


-fl 


0.11 


0.0 




1.09 


- 


-1 


0.19 


0.0 



The correlation energies are very similar for both 
modes and both nuclei with however a slightly larger gain 
for the Q32 mode. This small difference induces also an 
increase of dynamical deformations with respect to static 
ones. More significant are the differences obtained for 
the negative parity states whose excitations are lower by 
300 keV. Qualitatively similar predictions are obtained 
with the SLy4 interaction. 

In order to check the sensitivity of the GCM results 



on the magnitude of pairing correlations, we have per- 
formed calculations with pairing strengths that produce 
gaps larger or smaller by a factor 2. Qualitatively, the 
GCM results are not affected: the changes in the col- 
lective wave functions and dynamical deformations are 
marginal. The correlation energies vary by 10 — 20%. A 
doubling of the pairing gap results in approximately twice 
larger excitation energy for the negative parity states. 
We have found also that the results for ^*Zr are less sen- 
sitive to the pairing strength than for ^''Zr. 



IV. CONCLUSIONS 

We have studied two nuclei, ^°Zr and ^^Zr, which are 
doubly magic with respect to the tetrahedral symmetry, 
in order to establish whether they may serve as good 
candidates to obtain experimental evidence of a stable 
tetrahedral deformation. Our results, based on the mean- 
field approach, indicate that the spherical configuration 
of both nuclei is indeed unstable against octupole defor- 
mations. The presence of a tetrahedral minimum is the 
result of a delicate balance between shell effects and pair- 
ing. Due to the very shallow energy minima, quantum 
fiuctuations beyond mean-field play a significant role. In 
order to quantify this role, we have performed dynam- 
ical calculations using the GCM in two octupole direc- 
tions, specified by either tetrahedral or axial octupole 
deformations. The correlation energies associated with 
the octupole collective modes are large in both cases and 
lower significantly the energy of the spherical configura- 
tion. These results seem to be qualitatively independent 



of the Skyrme parametrization. 

We have performed the calculation in such conditions 
that each mode is clearly separated from the other, con- 
straining in particular the quadrupole deformations to 
be zero. Under these conditions it turns out that the 
correlation energies are slightly larger for the tetrahedral 
mode than for the axial one. Also the tetrahedral vibra- 
tion has smaller energy than the corresponding energy 
of axial octupole collective mode. Hence it was shown 
that thanks to dynamical effects, the tetrahedral mode 
is energetically more favorable as compared to the axial 
octupole mode. However the possible mixture between 
these modes cannot be excluded. 

In summary, our results suggest that both the axial 
and tetrahedral type of octupole correlations play an im- 
portant role in these nuclei, and they possess very similar 
characteristics. 
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